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Abstract 



A comprehensive, relativistic many-body approach to hadron structure 
is advanced based on the Coulomb gauge QCD Hamiltonian. Our method 
incorporates standard many-body techniques which render the approxima- 
tions amenable to systematic improvement. Using BCS variational methods, 
dynamic chiral symmetry breaking naturally emerges and both quarks and 
gluons acquire constituent masses. Gluonia are studied both in the valence 
and in the collective, random phase approximations. Using representative val- 
ues for the strong coupling constant and string tension, calculated quenched 
glueball masses are found to be in remarkable agreement with lattice gauge 
theory. 
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Our knowledge of the standard model cannot be considered complete until explicit gluonic 
degrees of freedom are found and understood ||T|. In an effort to address this issue we 
advance a comprehensive framework for consistently describing and understanding hadron 
structure - including the glueball and hybrid sectors. The model is motivated in part by our 
previous studies of relativistic and nonrelativistic || quark models and by the pioneering 
work of Le Yaouanc et al. ||, who, following the spirit of Nambu and Jona-Lasinio |§, 
have constructed a quark-based model of the QCD vacuum. Only a brief description of 
our approach is provided here. A full, detailed treatment will be contained in a separate 
communication. 

The idea is to build on the known successes of the constituent quark model for heavy 
quarks by considering a many-body relativistic Hamiltonian in a quasiparticle basis where 
dynamical chiral symmetry breaking and massive gluon modes are explicit. Such a model 
incorporates an extensive Fock space but reduces to the simple quark model in the valence 
approximation. Furthermore, the simultaneous presence of quark and gluon degrees of 
freedom permits studying their mixture in hybrid and glueball states. This is especially 
important since glueball searches tend to occur in meson-rich regions of the hadron spectrum 
and also because it may be years before lattice gauge calculations provide significant insight. 
This letter focuses on the gluonic sector of the model Hamiltonian, presenting the glueball 
spectrum calculation and a discussion of the associated approximation schemes. In the 
summary we comment on other issues regarding applications to mesons, baryons and hybrids. 

There have been a variety of previous glueball studies: the Bag Model QCD 
Sum Rules [PHI^I, the Constituent Glue Model Hl3| , and the Flux Tube Model fL4f . These 



approaches differ markedly in their mass predictions, sometimes by as much as 1 GeV, and 



no single approach has consistently reproduced lattice gauge measurements p5|-p^|. As 
stated above, our goal is to model QCD in a way which is in accord with the successes of the 
constituent quark model. Therefore we start from the Coulomb gauge QCD Hamiltonian H 
[HI and assume the existence of a set of phenomenological interactions H p h en such that H 
can be written as 



H = H + H I (1) 

where Hq is defined as 

H = K + Hphen, (2) 

and K stands for the kinetic energy 

K = |rfx^J(x)[-m-V + /?mJ^(x) + i|rfx[|E a (x)| 2 + |B a (x)| 2 ], (3) 

involving current quarks having masses m q and zero mass transverse gluons. The residual 
potential Hj = Hf CD - H phen is given by the difference between the original QCD and phe- 
nomenological interactions. The motivation for introducing phenomenological interactions 
is to generate a much weaker residual potential Hj at all energy scales. Furthermore, as 
in the phenomenological quark model, we will assume that in the quasiparticle quark and 
gluon basis this residual interaction can be approximated by the canonical QCD interaction 
Hamiltonian with a coupling g s that is small and saturates at low energies. Under this 
approximation, to any order in g s the residual interaction can in principle be derived using 
standard methods of time independent perturbation theory. However, since the infrared be- 
havior has already been determined from phenomenology, the residual interaction must be 
free from infrared divergences to avoid double counting. This can be achieved, for example, 
by imposing a cutoff, A IR , on Hj that removes coupling between quasiparticle states whose 
energy difference is smaller than the cutoff scale. The perturbative expansion of Hj then 
generates an effective Hamiltonian which has nonvanishing matrix elements below the cutoff 
and can either be added to H and diagonalized nonperturbatively or, because of the small 



coupling, included in the bound state perturbation theory using the eigenstates of H |L8|. 

For the applications considered in this work the phenomenological Hamiltonian is taken 
to be 

H phen = - l -j dxdyp a {*)V L (\* - y|)p a (y), (4) 

with color charge density p a (x) = *J(x)T a # 9 (x) + / abc A 6 (x)E c (x); V L is a linear confining; 
potential, 



Vb(|x - y|) = J^l x - y| (i - e- A ^" |x - y| ) (5) 

and N c = 3 is the number of colors. The string tension will be fixed at 0.18 GeV 2 , com- 
mensurate with Regge phenomenology and the naive quark model. Since this interaction is 
meant to represent soft physics, we have introduced an ultraviolet cutoff directly into the 
potential. Note that the approach advocated here requires that A p h en ~ Am and hence to 
completely describe low energy phenomena we must add the QCD interactions below the 
cutoff Aphen to this phenomenological Hamiltonian. To order a s we only keep the Coulomb 
potential and ignore self-energies, hyperfine interactions and vacuum corrections. This will 
be discussed in more detail below. The complete interaction in the Hamiltonian H, which 
we diagonalize nonperturbatively is thus given by Eq. (j|) with Vl replaced by 

V L (r) ->V(r) = V L (r) + V c (r), (6) 

where 

*fc(|x-y|) = -|^ ( 7 ) 

is the Coulomb potential. Finally we note that both Vc and Hi need to be ultraviolet 
regulated. The standard perturbative renormalization procedure may then be followed for 
these terms |[22|| . 

Regarding the structure of H v h en we have assumed that the bulk of the low energy 
dynamics of the qq and gluon-gluon systems may be described by a two-body interaction 
as shown and have employed a linear confinement potential in keeping with the constituent 
quark model and lattice gauge theory. Also note that the confining interaction is between 
color densities rather than scalar currents as is usually assumed in the constituent quark 
model. This is discussed below. Note that Eq.([|) implies that gluons may be confined into 
gluon-gluon bound states which forms the basis of our glueball investigation. However lattice 
gauge results |2(J indicate that static color octets become screened at very large distances 
and Eq.(^j) does not reflect this. Physically, one expects that the screening of the gluonic 
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confinement potential is due to the creation of low-lying glueballs at large valence gluon 
separation. Thus the model as presented here is similar in spirit to the naive quark model, 
where linear confinement in the qq sector is absolute and it is mixing to other Fock sectors 
which is responsible for the screening. Furthermore, we note that constituent gluons are not 
static, low lying glueballs tend to be compact, and that lattice gauge calculations support 
the picture of a diffuse spectrum. Finally, one should not think of the confinement potential 
as a flux tube because this quickly leads to conundrums about double counting gluonic 
degrees of freedom. The Coulomb gauge Hamiltonian makes a clear distinction between the 
gluonic color density and the Coulomb interaction - one which applies equally well to the 
confinement interaction. 

As in Ref. || the Fock space in which we diagonalize H is constructed from the vari- 
ational, BCS vacuum \Q) by an application of the constituent quark (antiquark) (D>) 
and gluon a^ creation operators with 

r dk r n 
= E J 7^1 A)S(k, A) + V(-k, A) J Dt(-k, A)] e *~ 



(27T) 

dk 



x 

A\r 1 

ikx 



a a (k) + a a (-k) t e ikx (8) 



W V^dki) 

and with the vacuum defined by B\Q) = D\Q) = a\Q) = 0. We note that the gluon operators 
are transverse so that one has 

[<(k), a}(q)t] = 5 ab (27r) 3 5(k - q) - kk) (9) 

The extra term in the final factor complicates the calculations of the glueball spectrum but 
is crucial to maintaining the correct gluonic degrees of freedom. 

The variational gap functions U, V and uo are obtained by minimizing the ground-state 
(vacuum) energy E = (fl\H \Q). This leads to gap equations which are equivalent to the 
Schwinger-Dyson equations for the quark or gluon self-energies in the rainbow approxima- 
tion. For the gluon spectral function uj the gap equation is given by 

„ (9) . s a(W =q * + ^J ^K (k + q) (1 + (fi . -"' k "^"""' (10) 



A similar variational treatment of the Hamiltonian in the quark sector results in the well 
known realization of dynamical chiral symmetry breaking by the BCS vacuum. 

The presence of the Coulomb term in V introduces a quadratic cutoff dependence, which 
can be removed by including the neglected terms in Hj (in particular the self-energy cor- 
rections resulting from the expansion of the Faddeev-Popov determinant and the transverse 
gluon exchange calculated to order a s ). However, the net effect of the order a s terms in Hi 
is expected to be small and we simply ignore them when solving the gap equation. Further 
investigations are in progress. 

A good fit to the numerical solution of Eq. (|T0|) is obtained with 

u{k) = ^fk 2 + m 2 g e- k2 / K2 . (11) 
If one defines the gluon mass in terms of the effective mass as m g = m g (0) where m g (k) = 



u>(k) 2 — k 2 then the proceeding fit yields m g = 0.8 GeV (and k = 13 GeV). 
The gluon condensate may be simply calculated within the context of the pairing ansatz. 
The result is 

£ GT G U = ZLll r (12) 

7T ^ 7T JO U!{K) 

We have allowed for the possibility that a s runs above the cutoff although this is not crucial. 
The calculated condensate agrees well with the QCD sum rule value of 0.012 GeV 4 |ZT| and 
is only weakly sensitive to the cutoff above A phen ~ 4 GeV. 

Just as in conventional nuclear structure theory, our BCS many-body vacuum state 
can be systematically improved by utilizing the Tamm-Dancoff (TDA), random phase 
(RPA), or even more accurate approximations involving exact diagonalization in an ex- 
tensive multiparticle-hole model space. In the glueball case we have performed both TDA 
and RPA calculations, however, because of the large constituent gluon mass we expect the 
Tamm-Dancoff approximation to be a reasonable one. Indeed the ++ glueball mass is 
shifted by less than 2% in going to the random phase approximation. In the glueball rest 
frame, the TDA gluon-gluon bound states are given by 
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\j pc ) = I ^y 3 4 pc (p)4(pM- P y\n) (is) 

with the glueball wave function xff C satisfying 

E^sr^ - -f / j^n\± + qi) 't?^ V c < k - q '^ C(k) 

-(h) + ^) «r + t / <l> k +«'>^^^ c H (14) 

Here the ^ PC 's are determined from coupling the two transverse gluons labeled by the 
Cartesian indices i,j to a state with total angular momentum J, parity P and charge 
conjugation C. For example 



^ ++ (k,q) = (l + (k-q)^ (15) 
7g- + (k,q) = ^ > (16) 

with more complicated expressions for J > 2. The functions are normalization matrices 
which arise from mixing between different LS states induced by the transverse nature of 
the gluon. We use the Coulomb and linear contributions to V in Eq. (14). There are no 
Faddeev-Popov terms and transverse gluon exchange is treated as a perturbation (it is of 
order l/m 2 g ). We note that it is not possible to make a two-gluon J = 1 state as is consistent 
with Yang's theorem. Such spurious states exist in models with explicitly massive gluons. 

There is an interesting property associated with divergences in the bound state equation. 
The linear potential is infrared divergent; however, this potentially problematic divergence 
is cancelled by the self-energy term in the kinetic energy. This cancellation appears to be a 
property of a density- density interaction. For example, the cancellation does not occur for 
scalar quark currents (and indeed, a stable vacuum cannot be obtained with an interaction 
between scalar currents). Furthermore, the cancellation appears to only occur in bound 
state equations for color singlet objects. This has been observed previously in the context 
of the Bethe-Salpeter approximation to the qq bound state problem |22f . 

The spectrum which results from numerically solving Eq. (|TJ|) is presented in Fig. 1 
along with recent results from lattice gauge calculations. The agreement is remarkably 



good, especially when it is recalled that the model has been completely fixed from qq phe- 
nomenology. The spectrum corresponds to the ansatz for a running a s (k) from Ref. |23| . 
Employing a fixed value for the strong coupling, a s = 0.4, produced no significant difference. 
Furthermore, the spectrum is essentially independent of the cutoff for A p h en ~ 4 GeV. We 
conclude that the model captures the essential features of glueballs. To our knowledge, this 
is the only model of gluonia which successfully reproduces lattice data and therefore it may 
provide important insight into glueball structure. 
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Fig. 1. Lattice Gauge and Model Glueball Spectra. The black bars are the results of the model 
calculation. Lattice results are indicated by open [16|, dark pTfl, and light [ITJ] boxes. 



Future, more comprehensive, studies may dictate phenomenological potentials beyond 
the two-body form used here. For example, the phenomenological 3 Pq decay vertex cannot 
be obtained from the density- density confining potential. Reconciling the naive quark model 
phenomenology with the model presented here should prove very instructive. Another issue 
is that the dynamical breaking of chiral symmetry should lead to a massless pion solution. 
Recall, however, that the BCS vacuum is not a true eigenstate of the Hamiltonian. Therefore 



S 



when diagonalizing the Hamiltonian in a truncated Fock space a chiral pion solution may 
not necessarily appear. However, in the random phase approximation which builds upon 
the BCS vacuum, the Nambu-Goldstone realization of chiral symmetry is preserved. In 



particular the Gell-Mann-Oaks-Renner J3] relation 



/>£ = -2m q {qq) (17) 



follows from Thouless' theorem |25[ applied to the chiral charge operator, Q$ = 
J dx 1 i r '(x)75 1 I r (x) via the expression 

2£ \(n\Q 5 \Q) RPA \ 2 (E n - E ) RPA = (Q\[Q 5 , [Q 6 ,H]]\Q) (18) 

n 

Examining the interplay of these issues will be instructive, especially as it has bearing on the 
rather mysterious nature of the hyperfme splitting and the ultimate utility of the potential 
quark model. 

Future work will also focus on baryon and hybrid structure. Towards this end, we have 
performed initial, but preliminary, calculations in the quark sector finding m g ~ 180 MeV 
and (qq) ~ —(100 Mev) 3 in agreement with Refs. PJ2^]. While m g is in rough agreement with 



phenomenology the low condensate value may be due to truncation of H phen to two-body 
form. We plan to extend this to higher terms and also incorporate the above mentioned 
many-body treatments. Of special interest will be an ambitious multi-particle/multi-hole 
diagonalization which will involve higher quark Fock state components. In particular the 
importance of such states as \qqq{qq)) for the proton will directly address the role of sea 
quarks and hidden flavor. Related to this is the insight this model provides concerning the 
proton spin. 

In summary, we have presented a unified, comprehensive approach to hadron structure 
based on non-perturbative relativistic field theory and the QCD Hamiltonian. The model 
provides both the appealing physical insight associated with the phenomenologically suc- 
cessful quark model and a consistent unified framework for studying issues such as chiral 
pions and quark-glue mixing. With the advent of CEBAF a challenging opportunity is at 
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hand to confront new precision data and to thoroughly investigate a wide variety of issues 
in hadronic physics. 
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FIG. 1. Lattice gauge and model glueball spectra. 
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